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ABSTRACT 

Close binary systems of compact stars, due to the emission of gravitational radiation, 
may evolve into a phase in which the less massive star transfers mass to its companion. 
We describe mass transfer by using the model of Roche lobe overflow, in which mass 
is transferred through the first, or innermost, Lagrange point. Under conditions in 
which gravity is strong, the shapes of the equipotential surfaces and the Roche lobes 
are modified compared to the Newtonian case. We present calculations of the Roche 
lobe utilizing the second order post-Newtonian (2PN) approximation in the Arnowitt- 
Deser-Misner gauge. Heretofore, calculations of the Roche lobe geometry beyond the 
Newtonian case have not been available. Beginning from the general N-body Lagrangian 
derived by Damour and Schaffer, we develop the Lagrangian for a test particle in the 
vicinity of two massive compact objects. As an exact result for the transverse-traceless 
part of the Lagrangian is not available, we devise an approximation that is valid for 
regions close to the less massive star. We calculate the Roche lobe volumes, and provide 
a simple fitting formula for the effective Roche lobe radius analogous to that for the 
Newtonian case furnished by Eggleton. In contrast to the Newtonian case, in which the 
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effective Roche radius depends only upon the mass ratio q = milm 2 , in the 2PN case 
the effective Roche lobe radius also depends on the ratio z = 2(mi + m 2 )/a of the total 
mass and the orbital separation. 

Subject headings: relativity — binaries: close — stars; mass loss 


1. Introduction 

During the evolution of a close binary system involving compact stars, the stellar separation 
shrinks due to the emission of gravitational waves. In the event that the stars are not of equal mass, 
and the less massive star has a larger radius than its companion, mass transfer may ultimately occur. 
Gravity wave emission generally causes the mutual orbit to circularize (Peters 1964). For circular 
orbits, conservative mass transfer can be modelled as Roche lobe overflow under the assumption 
that the star is not signihcantly disrupted due to tidal interactions. The Roche lobe is the innermost 
gravitational plus centrifugal equipotential surface encompassing both stars. 

In the model Roche lobe overflow, the radius of the less massive star is compared to the effective 
radius of its Roche lobe. Once the two radii become equal, because the Roche lobe radius decreases 
due to orbital decay, the star fills its Roche lobe and mass transfer occurs through the hrst, or 
innermost, Lagrange point Li. Lying on the Roche lobe, Li is located between the two stars on 
the axis connecting their centers and is also a saddle point of the gravitational plus centrifugal 
potential between the two stars. Due to its saddle point nature, the first Lagrange point acts as a 
gravitational funnel through which mass transfer occurs. 

Values of the Roche lobe radii as a function of orbital separation and mass ratio q = ra\lm 2 ^ 
where mi refers to the lighter star, have been tabulated by Kopal (1959) for the Newtonian case. 
Paczyhski (1971) and Eggleton (1982) have given analytical hts. We use Eggleton’s functional form, 
which has the advantage of being a continuous function of g, as a template in our work. 

In this work, we carry out calculations of Roche lobes beyond the Newtonian case. We employ 
the Arnowitt-Deser-Misner (ADM) form of post-Newtonian expansion and use the corresponding 
Lagrangian at the second order (2PN) level wherein terms up to (M/r)^, where M = mi + m 2 
and r is the distance, are retained. The same procedure as used in the Newtonian case for finding 
the Roche lobes is utilized. Our strategy is to (i) construct the effective potential for the point 
particle in the vicinity of two stars (the 3-body problem) in the co-rotating frame; (ii) evaluate 
equipotential surfaces and calculate the corresponding effective Roche volume and radius for this 
potential; and (hi) provide new htting formulae as Eggleton did for applications involving mass 
transfer. 

The organization of this work is as follows. In §2, we calculate the effective potential for three 
bodies at the 2PN level. We establish the Lagrangian in §2.1. The transverse-traceless part of 
the Lagrangian is evaluated explictly in §2.2 through the introduction of an approximation valid 
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for regions near mi for test particles. In §3, we evaluate the Roche lobes and their effective radii 
as a function of q and a relativity parameter for this potential, and provide a simple analytical 
fit. In this section, we also show the impact of post-Newtonian corrections on the positions of the 
Lagrange points and on the position of the center of mass. Our conclusions are contained in §4. 



Fig. 1.— The notation used in the evaluation of the Roche lobes in the 2PN approximation. 
Stellar masses are denoted by mi and m 2 and the point-particle mass is taken to be mo. Vectors 
Ryl {A = 0,1, 2) denote positions of the three bodies with respect to the origin O, r is the position 
of a generic point P, and is the position of this point with respect to the mass mA (we show 
only ro). The vectors Rab indicate positions of the three bodies with respect to each other. 


2. The 2PN potential for 3 bodies 

Roche lobes are defined through the acceleration that a point-like particle feels in the frame 
that is co-rotating with the two massive objects. While velocities of all three bodies disappear in 
this frame, accelerations are vanishing only for the two massive objects with masses mi and m 2 . 
The effective potential causes acceleration on the third object, the point-particle with mass mo. 
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Starting from the N-body Lagrangian in ADM coordinates derived by Damour & Schafer (1985), 
we obtain the 3-body Lagrangian for the situation depicted in Figure 1. We adopt the convention 
in which we denote masses with uppercase Latin indices {A, B,...) and coordinates with lowercase 
Latin indices {i,j, ■ ■ .)• Also, we use units such that G = 1 and c = 1. We express all inertial-frame 
velocities in terms of the rotating-frame velocities and the remaining rotationally induced part: 

VA = V™* -hw X FA, (1) 


where oo is the angular frequency of the rotating frame. Setting v™* = 0 for A = { 0 , 1 , 2 }, and 
v™* = 0 for A = { 1 , 2 }, but keeping v™* 7 ^ 0 enables us to find the acceleration on the point-like 
particle (body 0 ) from 
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where we have denoted the coordinate of the point-like body in the co-rotating frame by We 
have 


6 = 0, ^ = 0, and 6/0, (3) 

since Vq'’' = 6 We can find the effective potential for the particle 0 by separating out the “kinetic” 
part of the Lagrangian that contains terms that are quadratic in ^ and by treating the remaining 
part of the Lagrangian as the effective potential that we have to determine. It is straighforward to 
verify that this approach yields the Euler-Lagrange equations for 6 - After setting = 0, we use 
the resulting potential in order to trace the equipotential surfaces that correspond to the Roche 
lobes. 


2.1. The 2PN Lagrangian 


The computation of the effective Roche radii requires the effective potential that acts on a point 
particle in the vicinity of two massive bodies. In order to improve upon the existing Newtonian 
results, we utilize results that were obtained by using the post-Newtonian approximation of general 
relativity. 


The Roche problem requires a three-body Lagrangian in the case in which one of the bodies 
is a point-like particle of infinitesimal mass. Such results were derived for the more general N-body 
case by Damour &: Schafer (1985) who retained terms up to the second order (2PN) in M/r in the 
Arnowitt-Deser-Misner (ADM) coordinate gauge. 


For completeness, we list the main results of Damour &: Schafer (1985) here. The Newtonian, 
or zeroth order, result is familiar: 


1 


1 


Ln = i;^^ava + ^ 


ruAmB 


( 4 ) 


A,B^A 


rAB 
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where r^B = |R-a_b|- The first order post-Newtonian correction is 
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whereas the second order post-Newtonian contribution is 
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The term Utt refers to the transverse-traceless part of the Lagrangian potential which requires 
special treatment and is explicitly evaluated in § 2 . 2 . 

We now specialize to the case of three bodies with masses uiq, mi, and m 2 . For the Newtonian 
and the first order post-Newtonian correction, we obtain 


and 
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Because the second order post-Newtonian correction L 4 is small, and we are interested in 
computing the equipotential surfaces for a test particle, we assume 


mo <C mi, m 2 . 

We decompose the second order correction in order to facilitate its presentation: 


(9) 


L 4 = ^4“) + 0 (m 2 ) 


( 10 ) 


and we drop terms of O (m^). The decomposition is evident by comparing equation ( 6 ) with the 
following: 
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2.2. The transverse—traceless part of the Lagrangian Utt 

The calculation of the 3-body ADM Lagrangian from Damour & Schafer (1985) is somewhat 
lenghty, but straightforward, except for the transverse-traceless part of the interaction potential 
Utt in equations (6) and (14). It is not known in general how to evaluate this term explicitly, 
except for the two-body case (see Damour &: Schafer (1985); Ohta et al. (1973, 1974)). In order to 
circumvent this problem, and since we are interested only in the vicinity of the star with mass mi, 
we assume 


J’oi '^12 , and roi <C ro 2 . (22) 

and expand equations (5) and (6) in terms of roi/ri 2 . As in the rest of the Lagrangian, we assume 
the body 0 to be a point-particle and therfore drop all terms that are of quadratic (or higher) power 
in mg. 

Utilizing these two physically motivated assumptions, we expand Utt as 
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where A, B,C,D = 0,1,2 and g^i = dg{x)/dx^ for any function g{x). 
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where i = 1,2,3 denotes spatial components. Here, as in Figure 1, 


rA = r-R^, R^ = Rb+Ras, and ta = rB + 'Rab- ( 26 ) 

In the preceding equations, we denote the position of the integration point with r, the position (i.e. 
“trajectory”) of the body A with R^, the vector between the integration point and the body A 
with r^, and the vector defined by two bodies A and B with Rab- 

Two out of three terms in equation (23) are already available or easily derived from previous 
calculations of the two-body Lagrangian. The result for is given in Damour & Schafer 

(1985); Ohta et al. (1973, 1974): 
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It is straightforward to obtain an approximate result for by using an expansion in terms 

of roilri 2 and keeping only terms with nonpositive powers of rQijri 2 - This approximation has a 
straightforward physical interpretation that ro 2 ~ ri 2 around Ri. A quick calculation yields 
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However, the remaining term turns out to require a much more laborious calculation. We start 
from the expression 
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Instead of performing the entire calculation, we present all the necessary ingredients and techniques 
for an (alas, still long!) example. As an illustration, we show details of the calculation for the 
logarithmic term in equation (29). The procedures for the remaining terms are identical and are 
omitted here for brevity. 
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where terms ^(roi/ria) have been dropped. 

Various combinations of unit vectors can be expressed through derivatives with respect to 
particles’ distances about the integration point P 
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where di = d/dR\. For combinations of several analogous relations can be derived to be; 
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where in equations (32) through (34) we have abbreviated multiple derivatives as 
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After expressing the unit vectors in equation (30) as derivatives and taking derivatives in front 
of the integral sign, all integrations are performed over integrands that are combinations of powers 
of ri and r 2 only. Such integrals can be evaluated by the use of the formula (see Damour Sz Schafer 
(1985)) 


I{a,/3) 


(Px 

dvr 


Q (3 

^Ars = 


^r(g^)r(aa)r(-!a|±a) 

4 r{-f)r(-f)r(!^±f±5) 


^AB 


(36) 


After performing integrations, differentiations with respect to body trajectories Ra have to be 
performed. Integrations yield results that depend on Ra as powers of = |R^ — and in 

order to calculate we need to perform up to five consecutive differentiations with respect 

to Ra'- 
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d d d d d 1 


-{N + 4) ( + 5^^n% + 6^^n% 


+ 6^^n% + 


AB 


(40) 


+ + 5^^5^^ + n\B 

+ 

+ { 5 ^^5'^^ + + 5 *' 5 ^"^) n^B 

+ (^5*^,5*^™ + n(4B 

-(IV + 6) ^S^^n’fg + + 5^^n% 

+(5J'^n*]’g + + 5^^n% 

+5^^niB + + S^'^nfB 


+ (lV + 6)(iV + 8)n^“ 


ijklm I 


/• 


(41) 


where u^b'^ = n\B''^AB ■ • • ''^AS- differentiations with respect to R\ and R 2 can 

be related through the identity 


d 


d 


ani N 
^^B ^AB 


.N 


dR\ r^B 


(42) 


We now assemble these results. Firstly, after rewriting the unit vectors as derivatives through 
equations (31) and (32), we perform the necessary integrations by employing equation (36): 
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Finally, the application of differentiations in equations (37) through (41) yields the result 
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(45) 


12 


In order perform the full calculation for we need to compute the remaining terms in 

equation (29). They can be calculated by applying the technique described above on the following 
expressions: 
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A straightforward, but somewhat lengthy, calculation yields the following approximate result 
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—^ (2 ln(ri +ro + roi) - 

+tt4- (^0 nil + ni nil - n\ 

zrQi \ 


ri +ro rg(no • npi) 

m ‘^rh 

nil - n{ nil 



• noi) \ 
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The final result for the transverse-traceless term is 
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(51) 


3 . Roche Lobes for 2 PN 

We are now ready to find the equipotential surfaces that correspond to the Roche lobes and 
to compute volumes that are contained inside the lobes. We employ equations (11) through (21), 
supplemented with equation (51). After setting the three bodies on quasicircular orbits around 
the center of mass, we obtain the potential that gives rise to forces that act on the point-particle 
of mass mp. In Figure 2, we show the potential for y = 0 (we position the three bodies on the 
x-y plane so z = 0 automatically) and the corresponding Roche lobes for mass ratios q = 0.1, 0.2, 
0.5, and 1.0, respectively. Cusps on equipotential lines correspond to the first Lagrange points Li. 
All distances are scaled by the stellar radial separation a = ri 2 . Note that we only show the lobe 
around the first star, since our approximation (rpi <C ^ 12 ) is valid only in this region. 

As in the Newtonian case, the volume within the equipotential surface (Roche volume) grows 
with q for fixed a. However, the potential and the equipotential surfaces acquire an additional 
dependence. Unlike for the Newtonian case, the total mass modifies the result: for low q, the 
Roche volumes become smaller as the total mass increases, whereas for q greater than about 0.7 
the volumes increase. As for coordinate positions, where we have eliminated the separation a, we 
can introduce a new dimensionless parameter 

2M 
a ’ 

involving the ratio of the total mass M and the separation a. This parameter also corresponds to 
the ratio of Schwarzschild radius for M and the separation distance a. 

Integration of volumes enclosed by equipotential surfaces is straightforward; we utilize a 
Newton-Cotes type algorithm to find the enclosed volumes. Our results are shown in Figure 3 
for q G [0,1], and z = 0, z = 0.2, and z = 0.4, respectively. 
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Fig. 2.— Roche lobes and the corresponding potentials for y = 0. Coordinates x and y are scaled 
by the stellar separation a and are shown for q = 0.1, 0.2, 0.5, and 1.0, respectively. Results shown 
are for values of z = 0, z = 0.2, and z = 0.4, respectively. 
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Fig. 3.— Effective Roche lobe radii scaled by the stellar separation a versus g for z = 0, 

z = 0.2, and z = 0.4, respectively. 


In Figure 3, we also show results for the fitting function. Following Eggleton (1982), we choose 
the parametrization in which the scaled effective Roche radius is 


rRoche/a = Q{q)C{q,z), 


(53) 


where 


Qiq) = 


Oigq 


2/3 


/?Q (?2/3 + ln(l + gl/3) 


(54) 
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is the fitting function previously given by Eggleton (1982) and 

C{q, z) = 1 + 2 (ac - Pc) (55) 

is the correction function that stems from post-Newtonian effects. Fitting parameters of Q{q) are 
identical to the ones obtained by Eggleton (1982) (uq = 0.49 and Pq = 0.57) and values of 

ac = 1.951, and Pc = 1.812 (56) 

refer to the fitting function for C{q,z). This functional form describes extremely well the depen¬ 
dence of q for q < 1.0. We note that the crossover of reduced versus enlarged Roche lobes with 
respect to the Newtonian {z = 0) case occurs at q = {Pclctc)^ ~ 0.69, in which case the Roche 
radius is virtually ^-independent. 


3.1. Lagrange points and the center of mass 

While performing the computation of the effective Roche volumes and radii, it is necessary to 
find the first Lagrange point Li. In addition, there are two more extrema of the potential along the 
x-axis that correspond to the second and third Lagrange points, L 2 and L 3 , respectively. Moreover, 
it is necessary to find the position of the center of mass of the system while setting the particles 
onto quasicircular orbits. We briefly outline below how the positions of these points depend on the 
mass ratio and how the results get modified compared to the Newtonian case. We consider only 
the first three out of five Lagrange points (Li, L 2 , and L3) as L 4 and L 5 become local extrema 
only for nonvanishing velocity 


3.1.1. Lagrange points 

In the Newtonian case, the positions of the Lagrange points can be fully described in terms of 
the mass ratio q = Mi/M 2 . As the mass ratio drops towards 0, the position of the first Lagrange 
point quickly shifts away from the center of mass and in the direction of the lighter star. The 
positions of Lagrange points L 2 and L 3 show a much weaker g-dependence as can be seen in Figure 
4, where the center of mass is positioned at x = 0. The second point L 2 shifts towards the lighter 
mass for q < 0.2 after it passes through the maximum of its distance from the center of mass at 
q ~ 0.22. The third Lagrange point is slowly, but steadily, pulled toward the center of mass as 
g^O. 

As we might expect from the Roche lobe analysis, effects of post-Newtonian corrections give 
rise to an additional (Mi -I- M 2 )/ri 2 = M/a dependence. In our model, this can be parametrized 
with the quantity z = 2M/a. In Figure 4, we show two additional sets of lines for 2 = 0.2 and 0.4 
(2 = 0 corresponds to the Newtonian case). For all values of g given in Figure 4, the positions of Li 
and L 2 are slightly pulled towards the lighter mass compared to the Newtonian case. The position 
of L 3 is, however, pulled towards the center of mass for g > 0.45 and away from it for g < 0.4. 


- 19 - 



Fig. 4.— The positions of the Lagrange points Li, L 2 , and L 3 as a function of the mass ratio 
q = Ml/M 2 . Results are for z = 0, 0.2, and 0.4, respectively. 


3.1.2. Position of the center of mass 

In order to find the Roche lobes, we set the two stars on quasicircular orbits and set the center 
of mass at the origin. The post-Newtonian approximation modifies the position of the center of 
mass compared to the Newtonian result since the conservation of total momentum results in a 
slight shift of the position from the Newtonian case. In this section, we show how important these 
corrections are and parametrize them in terms of 2 : = 2M/a. 

We denote the ratio of the distance between mi and the center of mass by ai. The separation 
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between the two stars is a, and their ratio is 


— = P , and — = 1 — P . 
a a 


(57) 


Setting the center of mass at the origin yields the condition (de Andrade et al. 2001) 


Pl‘^ 1712 pi^ , m2 


Pmi < 1 + 5- - j + -5 


2mi^ 2a 8mi^ 4o V mi^ m2 


Pl^ P2^ 


+ 7 


Pi • P2 
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m 2 {mi + m 2 ) 




2m2^ 2a 8m2^ 4a \ m 2 ^ 


^ 

-5-^-?r + 7 


mim2 

Pi • P2 


mim2 


+ 


mi (mi + m2) 


(58) 


which contains an implicit dependence on P, since for quasicircular motion we can write the veloc¬ 
ities of the two stars as 


ui^ = P'^a^uP , and V 2 ^ = (1 “ P)‘^a^w‘^ , 

and their product as 

vi • V 2 = —/3(1 — P)oPuP , 

where at the 2PN level (Blanchet & Iyer 2003) 
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My .21 5 


io — —3 1 H-(z^ - 3)-I-I —I (— - 0 ^ + ^ 

a'^y a \a/48 

with V = q/{\ + q)'^. Squares of the two momenta are 

1 7 

pi^ = mi^P^a'^oj'^ + -mi^auP (2a + 5m2) -|- -mi^au;^m2/9 , 

2 2 4 4 4)4 

P2 = m2 O W 


(59) 


(60) 


(61) 


(62) 


P^ — Am 2 ^a'^ui'^P^ + (^m 2 ^a‘^Lo‘^ + ^m 2 ^au;^ (2a -|- Qa^uP + 5mi)^ /3^ 
(2a + QaPuP + 5mi) + ^m 2 ^aci;^ (—2a — 2aPuP — 12mi)^ P 
— -^m 2 ^auP (—2a — 2aPuP — 12mi) . 

The product of the two momenta can be computed to be 

Pi • P 2 = mim 2 a^u)^P^ — 2mim2P^a^u:^ + ymim 2 auP {pSaPuP -|- 4a -|- 5mi -|- 5m2) 


(63) 


1 7 

-|--mim 2 aa;^ (—4a — 2a^uP + 2 m 2 — 12mi) P — -mim 2 


2 2 
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(64) 


and their fourth powers (up to 2PN order) are 
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4 4 

mi ui 


(65) 
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P2 


4 


4 4 

m2 V2 . 


( 66 ) 


Inserting these results into equation (58), the position of the center of mass can be obtained 
by solving the 5^^ order polynomial 

0 = oo “1“ oi/3 “1“ Oj20^ ~\~ “t“ 0.4/3^ + 05/3^ , (67) 

with coefficients a, given by 

oo = —16mi^m2 — 16m2^mi — Sm2a^0J^ + 49mim2^a^u;^ — 6m2a^oJ^ — 16m2a^ 

+ 8 m 2 mia — 28 m 2 a^oJ^mi + 120mi^ 777 - 20 ^ 0 ;^ + 20 mim 2 a^u;^ , (68) 

oi = 32m,2^mi + 16m,2a^ + 16mia^ + 24m2a^w^ + 30m2a®a;'^ + 80m2a^a;^mi 

— 90771177120^14;^ + 327771^7772 — 9877717772^0^0;^ — 167T72777iO — 2307771^777,20^0;^ , (69) 

02 = 517771^77720^0;^ + 146777177720^0;'^ — 2477720^07^ — 8177717772^0^0;^ 

— 7277720^07^7771 — 6077720^07^ , (70) 

O 3 = 487772 O^ 0 ;^ 777 ,i + 8o^O;^777l + 6077720 ^ 07 ^ + 1077717772 ^ 0 ^ 0 ;^ + 107771 ^ 77720 ^ 07 ^ 

+ 877720'^0;^ — 84777177720^07^ , (71) 

04 = —3077720^07'^ — 20777177720^07'^ , (72) 

O 5 = 677720 ^ 07 ^ + 8777177720 ^ 07 ^ + 60 ® 07^7771 . (73) 


The results obtained by numerical root-finding are shown in Figure 5, in which we show the 
relative change of (5 compared to the Newtonian case: 

^/5 _ P - Pn 

P P ’ 

where Pn = 7772 /M. The dependence on the total mass M and on the separation a is entirely 
through their ratio and is parametrized in Figure 5 through the parameter z. 

The post-Newtonian analysis can be expected to be reliable up to moderate values of z (z < 
0.4 — 0.5). We note that in this regiont post-Newtonian corrections to /3 are less than about 4%. 
Moreover, in all cases of interest {q < 1.0), corrections to P are smaller than 0.5%. As expected, 
deviations from the Newtonian case increase for close or massive configurations, whereas for large 
separations or low masses the results converge toward the Newtonian case. Whereas for the more 
massive star (q > 1.0) the distance to the center of mass increases with 2 compared to the Newtonian 
case, the shift of the less massive star {q < 1.0) does not have a monotonic dependence on q. As 
expected from the symmetry of the problem, corrections vanish for q = 1.0. 


4. Conclusion 

We have utilized the second order post-Newtonian approximation in the Arnowitt-Deser- 
Misner gauge to calculate Roche lobe volumes. These results are an improvement over the Newto- 
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Fig. 5.— The position of the center of mass as a function of the relativity parameter 2 shown in the 
form of the ratio A/3//3 for g = 0.1... 10.0. We note an increase of [3 compared to the Newtonian 
case for q > 1.0. For q > 1.0, f3 decreases compared to the Newtonian case. 


nian case in that post-Newtonian gravity introduces corrections in the case of moderately strong 
gravitational field. 

In the course of our calculations, we have derived an approximate three-body Lagrangian that 
is valid in the case when one of the bodies is a point particle. This calculation requires an evaluation 
of the transverse-traceless term Utt of the Lagrangian for which an exact result is not available. 
However, as shown in this work, utilization of an approximation valid in the vicinity of the less 
massive star enables this problem to be circumvented. 

Using these results, we calculated Roche lobes in the 2PN effective potential in the co-rotating 
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frame and computed effective Roche lobe radii that can be used to model mass transfer through 
Roche lobe overflow. In addition, we computed changes to the positions of the Lagrange points 
and to the center of mass due to post-Newtonian effects. We find that corrections to Newtonian 
results for Roche lobe radii can be as significant as 20-30% at low mass ratio q < 0.1. Whereas for 
q > 0.7 the Roche lobe radius increases (~ 15% for q = 1.0), for low q’s the Roche lobe is smaller 
than in the Newtonian case. We have provided our results in the form of a simple htting formula 
that depends on two physical parameters; the mass ratio and the ratio of the total mass and the 
separation. 

Research support of the U.S. Department of Energy under grant number DOE/DE-FG02- 
87ER-40317 is gratefully acknowledged. 
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